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Abstract 

We employ quantum molecular dynamic simulations to investigate the behavior of benzene under 
shock conditions. The principal Hugoniot derived from the equation of state is determined. We 
compare our firs-principles results with available experimental data and provide predictions of 
chemical reactions for shocked benzene. The decomposition of benzene is found under the pressure 
of 11 GPa. The nonmetal-metal transition, which is associated with the rapid C-H bond breaking 
and the formation of atomic and molecular hydrogen, occurs under the pressure around 50 GPa. 
Additionally, optical properties are also studied. 
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I. INTRODUCTION 



The nature of hydrocarbons under extreme conditions is of great scientific interest and 
has recently attracted extensive studies. Specially, due to the possible importance in the 
formation of organic materials in the early solar system chemical reactions of small 
hydrocarbon molecules, such as acetylene, ethylene, methane, and benzene have gained 
particular attention. A major constitute of the "ice" layer in Uranus and Neptune, where 
the pressure could reach 600 GPa and temperature 8000 K, is thought to be methane 
and the planetary magnetic fields are greatly influenced by the electrical conductivity of 
methane at these conditions. Another important source of planetary carbon is thought to 
benzene, which has been found in the atmosphere of Jupiter and in carbonaceous chondrites 
. Naturally, comprehensive studies on the characteristics of benzene are indispensable 
for understanding the various astrophysical phenomena at high temperatures and pressures. 

Researches on benzene properties under shock conditions start from the plate impact 



experiments 



61] . As pressure increases, benzene has been found to decompose at 13 GPa 



along its Hugoniot 
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ll| . At high static and dynamic temperatures and pressures. 
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benzene mainly converts into defective carbon nanoparticles and H2, with the possib 
ucts of small concentrations of polycyclic aromatic hydrocarbons, such as alkanes 
Such reactions have also been found to be similar to the detonation of explosive materials 
14 1 . Nonmetal-metal transition of benzene was experimentally reported by Nellis et al. 
15l | and was attributed to thermal activation of defective carbon nanoparticles under the 
pressure between 20 and 40 GPa; further increase in pressure does not have much affect 
on the conductivities. To date, although a number of explanatory and predictive results in 
some cases have already been provided by experimental studies, less information is available 
for shock-introduced chemical reactions, and many fundamental questions of benzene under 
extreme conditions are still yet to be investigated. 

side, a tight-binding model has been applied to study the 



Concerning the theoretica 
properties of liquid benzene ij]. The results agree well with the experiments in (_D, u) 
diagram (shock velocity verse particle velocity), but the calculated temperatures remain too 
high. The thermodynamic properties of benzene were taken into account by classical force 
field method [l^]. While this method showed good agreement with the low-pressure exper- 
imental measurements, however, it is unsuitable in determining Hugoniot curve at higher 
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pressures. Overall, due to the intrinsic approximations of these methods, the electronic 
structure, which is predominant in determining the nonmetal-metal transition and equation 
of state (EOS) for benzene under shock conditions, is out of consideration. Thus, the change 
in the electronic structure of benzene under extreme condition is still recommended to be 
presented with a full quantum-mechanical description. On the other side, quantum molec- 
ular dynamics (QMD), where electrons are fully quantum mechanically treated, has been 
proven particularly suitable for the study of the chemical reactions, such as dissociation, 
ionization, and recombination of molecules. Meanwhile, the electronic structure, thermo- 
dynamical and optical pro per ties of warm dense matter have been successfully investigated 



by QMD simulations [18|, ll9||. According to the aspects mentioned above, the electronic 
structure, EOS and Hugoniot from density functional theory (DFT) are highly needed for 
shocked benzene. 

In the present work, we apply QMD simulations to study the properties of benzene along 
the principal Hugoniot. The EOS and the pair correlation functions (PCF) are determined by 
QMD simulations. The Kubo-Greenwood formula is used as a starting point for calculating 
the dynamic conductivity cr{uj), from which the dc conductivity is determined. The dielectric 
function e{uj) and reflectivity are then extracted. This paper is organized as follows. The 
simulation details are briefly described in Sec. [Tll The PCF, which is used to study the 
dissociation of benzene, and Hugoniot curve are given in Sec. Illlj In Sec. IIVI Nonmetal- 
metal transition and optical properties are discussed. Finally, we close our paper with a 
summary of our main results. 

II. COMPUTATIONAL METHOD 

In this study, we perform simulations for benzene by employing the Vienna Ab-initio 
Simulation Package (VASP) plane-wave pseudopotential code, which is developed at the 



Technical University of Vienna 
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2l| . Electrons are fully quantum mechanically treated 



through plane- wave, finite-temperature DFT 22|,|23|], and the electronic states are populated 



according to the Fermi-Dirac distribution at temperature Tg. The exchange correlation func- 
tional is determined by generalized gradient approximation (GGA) with the parametrization 
of Perdew-Wang 91 [24]. The ion-electron interactions are represented by a projector aug- 



mented wave (PAW) pseudopotential 25|, l26|]. Atoms move classically according to the 
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forces, which originate from the interactions of ions and electrons. The system is calculated 
with the isokinetic ensemble (NVT). In all the simulations, the system is kept in local ther- 
modynamical equilibrium by setting the electron (Tg) and ion (Tj) temperatures equal. The 
ion temperature Tj is kept constant every time step by velocity scaling. 

For molecular dynamic simulations, only F point is employed to sample the Brillouin 



zone, while 4x4x4 Monkhorst-Pack [27j scheme k points are used for the electronic struc- 
ture calculations. The plane-wave cutoff energy is set to be 600.0 eV. 48 carbon and 48 
hydrogen atoms (8 benzene molecules) are calculated in a cubic cell at separate densities 
and temperatures. The selected densities range from 1.2 to 2.7 g/cm^ and temperatures 
between 500 K and 6000 K to present the principal Hugoniot. All the dynamic simulations 
are lasted for 4~6 ps, and the time step selected for the atomic motion is 2 fs. Then, the 
system is equilibrated 200 steps and the properties are calculated via the data from the final 
300 steps. 

III. SHOCK EQUATION OF STATE AND PAIR CORRELATION FUNCTION 

The accurate calculation of electrical and optical properties depends on a precise de- 
scription of the materials properties, such as EOS. A crucial measurement of EOS data of 



benzene under shock condition is the Hugoniot [28|, which can be derived from conservation 
of matter, momentum, and energy for an isolated system compressed by a pusher at a con- 
stant velocity. Rankine-Hugoniot equation describes the locus of points in {E, P, V^)-space 
satisfying the relation as follows: 

(^0 - ^i) + liVo - Vi){Po + Pi) = (1) 

where E is the internal energy, V is the volume, P is the pressure, and the subscripts and 1 
refer to the initial and shocked state, respectively. In the canonical ensemble, shock adiabat 
between the initial and final states is described by the principal Hugoniot, which includes 
the locus of states {E, P, V) satisfying Eq. ([T]). In our present calculations, the pressure 
P is evaluated using the forces provided by VASP. The internal energy consists of the total 
energy from finite-temperature DFT calculation and zero-point energy. The determination 
of Hugoniot points is described as follows. For a given Vi, a series of simulations have been 
executed for different temperatures T. Then Ei and Pi are fitted to a cubic function of 
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T. The principal Hugoniot points {Ei, Pi, Vi) are evaluated by solving Eq. ([T]). For the 
present system, the initial density po is 0.874 g/cm^, and the internal energy Eq= 72.83 
eV/molecule at the temperature T=298 K. The initial pressure can be neglected compared 
to the high pressure of shocked states. Principal Hugoniot points derived from Eq. ([T]) are 
listed in Tab. [H 

Agreement between our calculated Hugoniot and experimental data is revealed in Fig. [H 
The region, which lies in the density range between 1.7 and 1.8 g/cm^, should be noticed. 
This region is accompanied with a rapid increase in pressure, which is attributed to the 
incipient dissociation of benzene and the formation of diamond-like nanoparticles. The 
transition of atomic structure results in the change of physical properties. Diamond is the 
hardest material, and the compression in such system will consequently lead to the rapid 
increase in pressure. Experiments indicate that benzene decomposes at 13 GPa 



D the rapid 



ll[ |. while our QMD simulations provide the value of 11 GPa. We further examine the EOS 
data at the densities ranging from 2.1 to 2.4 g/cm^ as shown in Fig. O The pressure shows 
a systematic behavior in terms of the density and temperature, except for the region (p=2.3 
~ 2.4 g/cm^ T=3500 ~ 4000 K, and P^50 GPa). Due to the rapid C-H bond breaking and 
the formation of atomic and molecular hydrogen, this region is featured by {dP/dT)y <0. 

TABLE I: Principal Hugoniot points derived from DFT-MD simulations at a series of density (/o), 
pressure (P), and temperature (T). 



p {g/cTD?) 


P (GPa) 


r(K) 


1.2 


4.3 


598 


1.4 


7.9 


814 


1.6 


10.0 


1093 


1.7 


10.9 


1200 


1.8 


19.3 


1324 


1.9 


19.2 


1437 


2.1 


33.7 


2963 


2.2 


43.2 


3370 


2.3 


54.4 


3749 


2.4 


67.2 


4334 
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FIG. 1: Comparison between the previous experimental data of Refs. 
QMD results for the principal Hugoniot. 



and the present 



The structural change of benzene under shock condition can be clearly reflected by PCF, 
which represents the possibility of finding a particle at a distance r from a reference atom. 
Typical C-C bond length in benzene molecule is 1.40 A, and the bond length of C-H is 1.09 
A in our calculations. In the case of diamond and H2, C-C and H-H have the bond length 
of 1.50 and 0.75 A, respectively. PCF, atomic structure, and charge density distribution 
at four densities of benzene along the principal Hugoniot are shown in Fig. [21 At p=0.874 
g/cm^ and T=298 K, one can see from Fig. [3]^a) that the PCF is featured by two peaks 
at 1.09 and 1.40 A, which represent the typical bond length of C-H and C-C in benzene 
molecule, respectively. These two peaks begin to reduce in amplitude and get broaden due 
to the thermal excitation at higher pressure [as shown in Fig. [31(b)]. For the densities p <1.7 
g/cm^ and temperatures T <1200 K, benzene turns out to retain in its molecular shape. 
This fact could be established through analysis of the atomic configuration and typical vr 
bond suggested by charge density. When p >1.7 g/cm^ and T >1200 K, diamond-like 
nanoparticles with vacancies and H substitutions are formed in the system. At this stage, 
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FIG. 2: P-T curves for four densities of benzene. 

C-C bond length changes into 1.50 A [Fig. [31(c)], which is a typical value of diamond, and 
the TT bond converts into sp hybrid bond. Benzene decomposes under the present condition, 
and the reference pressure is 11 GPa. With further increase in pressure, the dissociation of 
C-H bond becomes remarkable, and the formation of H2 molecule is suggested by the peak 
at r=0.75 A, as shown in Fig. [3]^d). 

IV. DYNAMIC CONDUCTIVITY AND OPTICAL PROPERTIES 

The real part of the dynamic conductivity ai{uj) is derived from the Kubo-Greenwood 
formula: 



2^p2fc2 N N 3 



(2) 




FIG. 3: Calculated pair correlation function for C-C (black line), C-H (red line), H-H (green line) 
at four densities of benzene along the principal Hugoniot. The atomic structure, where carbon and 
hydrogen atoms are denoted by gray and white balls respectively, and the relative iso-surface of 
charge density (blue regimes) are also provided in the insets, (a) p=0.874 g/cm^, T=298 K; (b) 
P=1A g/cm^, r=814 K; (c) p=2.1 g/cm^, r=2963 K; (d) p=2.3 g/cm^, r=3749 K. 

where /(ej,k) describes the occupation of the ith band, with the corresponding energy ej^k 
and the wavefunction ^ at k, and w(k) is the k-point weighting factor. The imaginary 
part of the dynamic conductivity o'2{uj) follows from the Kramer-Kronig relationship. Then, 
the real and imaginary parts of the dielectric function are given by ei{uj) = 1 — ■^a2{uj) 
and £2(^7) = -^ai{uj), respectively, from which the optical constants such as the refractive 
index n{uj), extinction coefficient k{uj) and reflectivity r{uj) may now be derived (ei=n^— /c^, 
e2=2nk, and r={i^^|J±g). 

The behavior of o"i(u;) at different Hugoniot points is summarized in Fig. HI Overall, we 
flnd that the real part of dynamic conductivity [cri(co')] is featured by a main peak around 
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FIG. 4: ai{uj) along the principal Hugoniot. 



10.0 eV and approach to vanish for the energy higher than 25.0 eV. For the pressure lower 
than 30 ~ 40 GPa (p <2.1 ~ 2.2 g/cm^), the dc conductivity, which could be given as: adc = 
lim cri(co'), is zero, which obviously shows the insulating nature of benzene at low pressure 
regime. With the increase of pressure along the principal Hugoniot (p=2.3 ~ 2.4 g/cm^), 
three prominent new features occur in the frequency-dependent conductivity: (i) The main 
peak in ai{uj) moves towards lower energies; (ii) There tends to develop an additional little 
peak around the energy of 2.5 eV, which again signifies the occurrence of diamond-like 
structure in the system; (iii) The dc conductivity is no longer zero. This indicates that 
nonmetal-metal transition, which is associated with the rapid C-H bond breaking, happens at 
the high-pressure Hugoniot points around 50 GPa (p~2.3 g/cm^ and T ~3700 K), consistent 
with the predicted value derived from experiment measurement 15j- Furthermore, the 
optical reflectivity of benzene along the principal Hugoniot is shown in Fig. [5l The pressure- 
induced change in reflectivity can be clearly indicated in the energy range of 1.0~5.0 eV 
(the reference wave length is correspondingly from 250 to 1250 nm). For the photon energy 
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FIG. 5: Optical reflectivity of benzene along the principal Hugoniot. 

£'=1.0 eV, a considerable increase in reflectivity (from 0.2 to 0.6), which is attributed to 
high pressure induced nonmetal-metal transition, is observed along the principal Hugoniot. 



V. CONCLUSION 

In summary, the QMD simulations presented here demonstrate the properties of shock 
compressed benzene. The EOS and principal Hugoniot are determined according to the 
simulated results, and the results show excellent agreement with the available experiments. 
Our calculations have shown that the benzene molecule begins to dissociate around 11 GPa. 
The considerable C-H bond breaking appears at the pressure around 50 GPa. The nonmetal- 
metal transition also occurs at this stage, which will lead to the increase in optical reflectivity 
from 0.2 to 0.6. 



10 



Acknowledgments 



This work was supported by the Foundation for Development of Science and Technology 
of China Academy of Engineering Physics under Grant No. 2009B0301037. 



[1] M. L. Elert, S. V. Zybin, and C. T. White, J. Chem. Phys., 118, 9795 (2003). 

[2] W. B. Hubbard, Science, 214, 145 (1980). 

[3] D. J. Stevenson, Annu. Rev. Earth Planet Sci., 14, 257 (1982). 

[4] K. Mimura, M. Kato, and R. Sugisaki, Geophys. Res. Lett., 21, 2071 (1994). 

[5] K. Mimura, Geochim. Cosmochim. Acta, 59, 579 (1995). 

[6] R. Dick, J. Chem. Phys., 52, 6021 (1969). 

[7] G. E. Hauver, J. Appl. Phys., 36, 2113 (1965). 

[8] R. D. Dick, J. Chem. Phys., 52, 6021 (1970). 

[9] R. H. Warnes, J. Chem. Phys., 53 1088 (1970). 

[10] O. B. Yakusheva, V. V. Yakushev, and A. N. Drcmin, High Temp. -High Press., 3 261 (1971). 

[11] N. C. Holmes, G. Otani, P. McCandless, and F. H. Ree, High Press. Res., 5, 669 (1990). 

[12] W. J. Nellis, F. H. Ree, R. J. Trainor, A. C. Mitcheh, and M. B. Boslough, J. Chem. Phys., 
80 2789 (1984). 

[13] K. Mimura, M. Kato, and R. Sugisaki, Earth Planet. Sci. Lett., 133 265 (1995). 

[14] L. L. Davis and K. R. Brower, J. Phys. Chem., 100 18775 (1996). 

[15] W. J. Nellis, D. C. Hamilton, and A. C. Mitcheh, J. Phys. Chem., 115 1015 (2001). 

[16] S. R. Bickham, J. D. Kress, and L. A. Collins, J. Phys. Chem., 112 9695 (2000). 

[17] J. B. Maillet and N. Pineau, J. Phys. Chem., 128 224502 (2008). 

[18] T. R. Mattsson and M. P. Dcsjarlais, Phys. Rev. Lett., 97 017801 (2006). 

[19] A. Kietzmann, B. Hoist, R. Redmer, M. P. Desjarlais, and T. R. Mattsson, Phys. Rev. Lett., 

98 190602 (2007). 

[20] G. Kresse and J. Hafner, Phys. Rev. B, 47 R558 (1993). 

[21] G. Kresse and J. Furthmiiller, Phys. Rev. B, 54 11 169 (1996). 

[22] T. Lenosky S. Bickham, J. Kress, and L. Collins, Phys. Rev. B, 61 1 (2000). 

[23] S. Bagnier, P. Blottiau, and J. Clerouin, Phys. Rev. E, 63 015 301(R) (2001). 



11 



[24] J. P. Perdew, Electronic Structure of Solids (Akademie Verlag, Berlin, 1991). 

[25] P. E. Blochl, Phys. Rev. B, 50 17953 (1994). 

[26] G. Kresse and D. Joubert, Phys. Rev. B, 59 1758 (1999). 

[27] H. J. Monkhorst and J. D. Pack, Phys. Rev. B, 13 5188 (1976). 

[28] Y. Zeldovich and Y. Raizer, Physics of Shock Wave and High Temperature Hydrodynamic 
Phenomena (Academic Press, New York, 1966). 



12 



